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energies from 300 MeV to 3 GeV. Differential cross sections da/d£2, da/dt, and the total cross section are 
presented. For < 3 GeV, the angular distributions agree well with the SAID parametrization. At 
photon energies above 1.5 GeV, a strong forward peaking indicates t-channel exchange to be the dominant 
process. The rapid variations of the cross section with energy and angle indicate production of resonances. 
An interpretation of the data within the Bonn-Gatchina partial wave analysis is briefly discussed. 

PACS: 13. 30. -a Decays of baryons, 13.60.Le Meson production 14-20. Gk Baryon resonances with S = 



Received: February 5, 2008/ Revised version: 
Abstract. Photoproduction of 7r° mesons was studied with the Crystal-Barrel detector at ELSA for incident 



1 Introduction 

Due to their substructure, nuclcons exhibit a rich spec- 
trum of excited states. A survey of the resonances ob- 
served so far can be found in [1]. In spite of considerable 
theoretical achievements, attempts to model the nuclcon 
spectrum with three constituent quarks and their interac- 
tions still fail to reproduce experimental findings in im- 
portant details. In most quark-model based calculations, 
more resonances are found than have been observed ex- 
perimentally [2,3]. However, the quark model is only an 



approximation and may ovcrpredict the number of states 
[4,5]. Alternatively, these missing resonances could as well 
have escaped experimental observation due to a weak cou- 
pling to N7r which makes them unobservablc in clastic 7rN 
scattering. 

Resonances with small N-7T couplings are predicted to 
have sizable photocouplings [2]. Thus, photoproduction of 
baryon resonances provides an alternative tool to study 
nucleon states. New facilities such as ELSA 1 at Bonn, 
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Graal 2 at Grenoble, Jefferson Lab (Virginia), MAMI C 3 
at Mainz, and SPring-8 (Hyogo) offer the opportunity to 
investigate photoproduction for E~ t > 1 GeV and to study 
nuclcon resonances above the first and second resonance 
region. 

Good angular coverage is needed to be able to extract 
the resonant and non-resonant contributions in a partial 
wave analysis. The Crystal-Barrel detector at ELSA is 
thus an ideal tool for studying nucleon resonances. 

Here, we present differential cross sections for the re- 
action 

7P^P7T°. (1) 



Nucleon resonances contributing to this reaction are not 
expected to belong to the class of missing resonances but 
when searching for these, identification of known reso- 
nances in photoproduction is an important step. 

The results on ir° and r\ photoproduction using the 
CB-ELSA detector were communicated in two letters [6, 
7]. In this paper, we give full account of the experiment 
and of data analysis of the reaction 7p — ► p7r°. A publica- 
tion covering all aspects related to r\ photoproduction is 
in preparation [8]. 

This paper is organized as follows: In section 2, we 
give a survey on the data already published before the 
CB-ELSA experiment. The experiment itself is then de- 
scribed in section 3. Section 4 provides a detailed descrip- 
tion of the data taken during the first data taking peri- 
ods and of the methods used in the event reconstruction. 
The determination of the differential cross sections and 
the treatment of systematic errors are discussed in sec- 
tion 5. Section 6 contains a short description of the PWA 
method used to extract the contributing resonances from 
the data [9—11]. A summary is given in section 7. 



2 Previous Results on 7r° Photoproduction 

First data on photoproduction of neutral pions date back 
to the 1960's. The older experiments were limited in an- 
gular coverage and energies. Fig. 1 shows the available 
differential cross section data for the energy range from 
0.3 to 3 GeV. The data are from the SAID 4 compilation 
[12]. 

In general, there is good coverage at lower energies. At 
beam energies above 2 GeV, there are only few data points 
measured by different experiments at discrete angles. 

The description of the data using the SAID model is 
also shown in Fig. 1. The SAID model was fitted to these 
data points and thus, is only reliable up to an incoming 
photon energy of about 2 GeV. 

In the following, a summary of the experiments per- 
formed and published after 1979 is given. 

Yoshioka et al. [13] extracted differential cross sec- 
tions for photon-beam energies between 390 and 975 McV 
in energy bins of 20 to 25MeV. They measured at 11 



2 Grenoble Anneau Accelerateur Laser 

3 MAinz MIcrotron 

4 Scattering Analysis Interactive Dial-In Program 



do/dQ |>b/sr] 




1.5 
1 

0.5 


1.5 

1 

0.5 



1300- 1350 



1550-1600 



1800- 1850 



2050-2100 



2300 - 2400 



1600 - 1650 




1850 - 1900 



2100-2150 



2400 - 2500 



1650-1700 



1900-1950 



2150-2200 



2500 - 2600 



1700 - 1750 



1950-2000 



2200 - 2250 



2600 - 2800 



1750 - 1800 



2000 - 2050 



2250 - 2300 



2800 - 3000 



I ' T ' I ' T 1 I 1 'I' 1 I 1 T ' I ' T ' I ' 'I' ' I ' T ' I ' 'I 
-1 -0.5 0.5 1 -0.5 0.5 1 -0.5 0.5 1 -0.5 0.5 1 -0.5 0.5 1 



cos0 cm 

Fig. 1. Differential cross sections da/dQ versus cos<9 C m as 
measured by various experiments before CB-ELSA. The data 
stem from [12]. The SAID model for the respective energy in- 
terval is shown as solid line. The photon energy range in MeV 
is indicated in the subfigures. 
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scattering angles covering a range in the center-of-mass 
system (cms) of 15° < cnl < 130° corresponding to 
-0.64 < cos(9 cm < 0.97. 

At the Saskatchewan Accelerator Laboratory, differ- 
ential cross sections were measured for 11 energies within 
25MeV above the 7r°-photoproduction threshold [14]. The 
Igloo spectrometer was used, which was especially de- 
signed for an excellent 7r°-detection efficiency. A full an- 
gular coverage from 0° to 180° was achieved. 

Beck et al. measured differential cross sections at the 
electron accelerator MAMI for five energy bins from thresh- 
old at 144 MeV up to photon-beam energies of 157 McV 
[15] and for six further bins between 270 MeV and 420 MeV 
[16]. Both experiments used a linearly-polarized photon 
beam produced via coherent brcmsstrahlung. In [16], the 
reaction was studied with the DAPHNE 5 detector which 
covers ~ 94 % of the solid angle. 

Differential and total cross sections were determined 
with TAPS 6 at MAMI [17]. The TAPS collaboration mea- 
sured at nine different energies between threshold and 
280 MeV. The setup used five blocks of crystals resulting 
in coverage of the full angular range in <9 cm , but only in 
partial coverage of the azimuthal angle. The experiment 
covered about 50% of the total solid angle. The published 
data are restricted to energies below 152.5 MeV. The con- 
tribution of the Eq + multipole was extracted and com- 
pared to predictions from chiral perturbation theory and 
low-energy theorems. 

Krusche et al. [18] extended the experiment to cover 
the energy range up to 792 MeV, which was the maxi- 
mum possible energy at MAMI. A comparison was made 
between ir° photoproduction off protons and off deuterons. 

Schmidt et al. [19] measured differential and total cross 
sections for ir° photoproduction for incoming photon ener- 
gies between threshold and 165 MeV. The measured data 
points were then also compared to predictions based on 
chiral perturbation theory and low-energy theorems. 

Ahrens et al. [20] measured differential cross sections 
for 12 energies in the photon-energy range between 550 
and 790 MeV using the DAPHNE detector at MAMI. The 
availability of a polarized photon beam and a polarized 
target made it possible to measure the helicity difference 
(T3/2 — 0"i/2 ■ The goal of the experiment was to test the 
GDH 7 sum rule [21]. 

At BNL 8 , photoproduction cross sections for tt° mesons 
were measured using LEGS 9 . Final-state particles were de- 
tected in an array of six Nal crystals. The data cover the 
beam-energy range from 213 MeV to 333 MeV. Unpolar- 
ized differential cross sections as well as beam asymmetries 
were determined [22]. 

This compilation shows that (almost) all data taken af- 
ter 1979 cover only the lower photon-beam energy range 
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up to 1 GeV. Data above 1 GeV stem from even older ex- 
periments; they are included in the SAID database [12], 
too. The only recent publications comprising new exper- 
imental data on ir° photoproduction originate from the 
GRAAL Collaboration [23] and from this experiment [6]. 

All available data (except [6,23]) were reproduced well 
by the SAID model solution SM02 [12]. Data and fit re- 
sults are shown in Fig. 1. The partial- wave analysis (PWA) 
also included 7rN-scattering data and determined masses, 
widths, and photocouplings of baryon resonances. The fit 
covered a photon-energy range up to 2 GeV. Contributions 
of the following resonances were extracted [24] : 

N(1440)Pn, N(1520)Di 3 , N(1535)Sn, N(1650)Sn, 
N(1675)Di 5) N(1680)Fi 5) Z\(1232)P 33 , Z\(1620)S 3 i, 
Z\(1700)D 33 , Z\(1905)F 35 , Z\(1930) D 35 , and Z\(1950)F 37 . 

Recently, a new SAID model solution (SM05) has been 
released which, in addition, takes into account the data 
presented in this paper, our data on 7P — > pry [7] and the 
new GRAAL data on 7p — ► p7r° [23]. 

The Mainz unitary isobar model MAID concentrates 
on the photon-beam energy range below 1 GeV [25]. Ex- 
tensions of MAID are in preparation [26]. Besides SAID 
and MAID there is a large number of other approaches to 
describe photon- and pion-induced production of mesons; 
a survey is given in the recent paper by Matsuyama, Sato 
and Lee [27]. 



3 Experimental Setup 

3.1 Overview 

The CB-ELSA experiment was performed at the electron 
stretcher accelerator ELSA. The maximum achievable elec- 
tron energy is 3.5 GeV. 

The accelerator complex comprises three stages. A lin- 
ear accelerator (LINAC) pre-accelerates electrons emit- 
ted by a thermionic electron gun to an energy of 20 MeV. 
These are transferred to a booster synchrotron, where they 
can be accelerated up to 1.6 GeV forming a pulsed beam 
(/ = 50 Hz). The electron bunches are injected into the 
stretcher ring over several cycles. The stretcher ring accu- 
mulates these bunches and accelerates the electrons up to 
the required energy. The electrons are made available for 
the CB-ELSA experiment via slow resonance extraction, 
with a duty factor of nearly 90%. 

The experimental setup is shown in Fig. 2. A first short 
overview of the different detector components is given be- 
low, which are then discussed in more detail in the follow- 
ing subsections. 

The electrons hit a radiator target, where they pro- 
duce bremsstrahlung. The energy of the photons, in the 
range between 25 % and 95 % of the primary electron en- 
ergy, was determined via the detection of the correspond- 
ing scattered electrons in the tagging system. The primary 
electron beam of unscattered electrons was stopped in a 
beam dump situated upstream of the Crystal Barrel. The 
photon beam hit the liquid hydrogen (LH2) target (length: 
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Fig. 2. Experimental setup CB-ELSA, Bonn. The electron beam enters from the left side, hits the radiator, and produces 
bremsstrahlung. The photons are energy tagged and hit an LH2 target in the center of the Crystal Barrel. Charged particles 
leaving the target are identified in the inner scintillating-fiber detector, photons are detected in the CsI(Tl) calorimeter. A 
photon counter for the flux determination is placed further downstream and is not shown in the figure. 



I = 52.84 mm, diameter: d = 30 mm) placed in the center 
of the Crystal-Barrel detector. 

The Crystal-Barrel detector forms the central compo- 
nent of the experiment. It consists of 1380 CsI(Tl) crys- 
tals and has an excellent photon-detection efficiency. The 
large solid-angle coverage and the high granularity allow 
for the reconstruction of multi-photon final states. A more 
detailed description of the Crystal-Barrel detector can be 
found in [28]. 

Charged particles leaving the target cannot be un- 
ambiguously identified by their energy depositions in the 
calorimeter. They were detected in a three-layer scintillat- 
ing fiber detector surrounding the target [29]. The first- 
level trigger of the experiment exploited detection of pro- 
tons; it was provided by the tagging system in coincidence 
with the fiber detector. For the photon-senstitive second- 
level trigger a FAst Cluster Encoder (FACE), based on cel- 
lular logic, provided the number of clusters in the Crystal 
Barrel. A segmented total-absorption oil Cerenkov counter 
(not shown in Fig. 2) was placed further downstream to 
determine the total photon flux traversing the target. 

3.2 Tagging System and 7-Counter 

The photon beam was produced by bremsstrahlung off an 
amorphous copper foil of 3/1000 radiation length thick- 
ness. The total rate in the tagging system during the 
beamtime was 1-3T0 6 Hz. The unscattered electron beam, 
deflected by 7°, was annihilated in a beam dump con- 
sisting of heterogeneous materials including lead, boron- 
carbide, and polyethylene. 

The tagging system is shown in Fig. 3. It serves to anal- 
yse the electron momentum spectrum behind the radia- 
tor. According to their energy loss, electrons are deflected 
in the dipole field of the tagging spectrometer; the field 
strength (1 - 2T) was set according to the chosen beam 
energy. Assuming single photon radiation, the photon en- 
ergy Ery is given by the beam energy E e ^ elsa and the 



energy of the detected electron E' e , 

E~f = E Ct elsa — E' c . (2) 

Tagged photons arc assigned to hadronic events in the 
Crystal-Barrel setup by coincidences. Thus, the tagging 
system consists of two distinct parts, a 4 cm thick scin- 
tillator array to provide fast timing information, and two 
MWPCs 10 with a total of 352 channels (with four chan- 
nels overlap). Their spatial resolution translates into an 
energy interval per channel of 0.1 (0.5) MeV at the high- 
est E 1 and 10 (30) MeV at the lowest E 1 for the beam 
energies of 1.4 (3.2) GeV, respectively, which were used 
for our measurements. The scintillator array was read out 

10 Multi-Wire Proportional Chambers 



Two MWPCs 




Fig. 3. Layout of the tagging detector: The ladder of 14 scin- 
tillation bars and two proportional wire chambers are shown. 
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via photomultiplicrs. A logical OR of the left-right coinci- 
dences from all scintillators was required in the first-level 
trigger. 

The tagging system was calibrated by direct injection 
of a very low intensity e~-bcam of 600 MeV or 800 MeV, 
after removing the radiator. Variation of the magnetic field 
of the tagging dipolc enabled a scan of several spatial posi- 
tions over the MWPCs. For a given wire, electron momen- 
tum and magnetic field strength are proportional. As long 
as saturation effects can be ignored, the magnetic current 
is proportional to the field. 

The calibration was checked by Monte-Carlo trace sim- 
ulations of the electron trajectories through the tagging 
magnet. Geometry of the setup, dimensions of the electron 
beam, angular divergences, multiple scattering and M0ller 
scattering in the radiator foil and the air were taken into 
account. From these simulations, the energy of each of 
the MWPC wires was obtained by a polynomial fit. The 
uncertainty of the simulation was estimated to be of the 
same order of magnitude as the energy width of the re- 
spective wire. Deviations between the 600MeV/800MeV 
calibration could be attributed to effects of magnetic field 
saturation. Hence, the calibration of the highest photon 
energy points relied on an extrapolation of the simulated 
trajectories. 

The hit-wire distribution of the MWPCs was measured 
with a minimum-bias trigger at a fixed rate of 1 Hz. This 
trigger required only a hit in the tagging system and was 
thus independent of hadronic cross sections. Accepted wire 
hits had to be isolated (one hit or cluster of hits); the small 
background was estimated from the time distribution of 
the associated tagger scintillator. 

The absolute normalization was treated as a free pa- 
rameter which was determined by fitting the measured 
angular distributions to the SAID cross sections. For the 
low-energy data (-E e , elsa = 1.4 GeV), the normalization 
constant was determined for each energy bin, and an error 
of ±5% was assigned to it. The E c , elsa = 3.2-GcV data 
cover a region for which no SAID prediction is available. 
An energy-independent normalization constant was deter- 
mined from a comparison of our differential cross sections 
with SAID for 0.8 < E 1 < 1.7 GeV. A systematic error 
of ±15% is estimated to account for possible variations of 
the background across the tagger. 

Downstream of the tagger and behind the Crystal- 
Barrel detector a total absorption oil Cerenkov photon 
counter was mounted. It consists of three segments, each 
made of a hollow lead cylinder with 10 lead blades sur- 
rounded by mineral oil. The light produced by traversing 
particles was detected by two photomultipliers per seg- 
ment. 



3.3 Liquid H 2 Target 

Fig. 4 shows a side view of the assembly, the gas liquefier, 
one half of the Crystal-Barrel calorimeter, target cell, and 
inner scintillation fiber detector. The target cell cylinder 
is made of 125 /im Kapton foil; entrance and exit window 
have a thickness of 80 /xm. The massive cold head of the 



1.1 m 




1.4 m 

Fig. 4. The liquid hydrogen target with target cell (1), inner 
SciFi detector (2) and photomultiplier readout (3) are situated 
in the center of the Crystal-Barrel calorimeter (4). Due to the 
geometry, there is a large distance between liquefier (5) and 
target. 



H2 refrigerator was mounted 2.5 m outside of the detector 
to avoid obstruction of the detector acceptance. 

The Air Products (model CSA-208-L) refrigerator con- 
sists of a two-stage cooler head operating in a closed he- 
lium circuit according to the Gifford-McMahon principle. 
In the first stage, H2 is liquefied. The second circuit con- 
sists of a large H2 cold gas reservoir and the target cell 
connected via two Kapton pipes of 75 /xm foil thickness 
serving as liquid H2 supply and gaseous H2 removal pipes. 



3.4 Crystal Barrel 

The Crystal-Barrel calorimeter is designed to provide high- 
efficiency photon detection with good energy and spatial 
resolution over an energy range from 20 MeV to 2 GeV. 
All crystal modules are arranged in a vertex-pointing ge- 
ometry forming the shape of a barrel. The Crystal Barrel 
consists of 1380 CsI(Tl) crystals with a length of 30 cm 
corresponding to 16 radiation lengths and covering almost 
the complete solid angle. Each of the crystals is wrapped 
with a 0.1mm titanium foil for mechanical stability and 
protection and a 2 mm support foil at the end (see Fig 5). 

Fig. 6 shows a schematic diagram of the crystal ar- 
rangement. The special shape of the calorimeter requires 
13 different types of crystals. They are arranged in 26 rings. 
Each ring has either 60 or 30 crystals, i.e. a single crystal 
covers 6° (rings 1-10) or 12° (rings 11-13) in azimuthal an- 
gle <P. In polar angle, a range from 12° to 168° is covered 
corresponding to a solid angle of 97.8% of 47r. 
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Fig. 5. A crystal module: titanium case (1), wavelength shifter 
(2), photo-diode (3), preamplifier (4), optic fiber (5), case cover 
(6). 
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Fig. 6. Mounting scheme of the Crystal-Barrel calorimeter. 
The numbers indicate the different crystal shapes. 



The size of the crystal segments was adjusted to the 
Moliere radius of Csl crystals, ru = 3.5 cm. The 12° open- 
ing on either side of the barrel is necessary for technical 
reasons. The segmentation limits the spatial resolution to 
20mrad (« 1.1°) in $ and 0. However, it allows the sepa- 
ration of two photons stemming from the decay of a ir° me- 
son with a maximum momentum of 1 GeV/c correspond- 
ing to a minimum opening angle of 16.6°. The energy res- 
olution of the calorimeter is empirically described by 



E 



2.5% 



</E [GcV] 



(3) 



For the readout of the scintillation light, silicon photodi- 
odes are used. They are placed on a wavelength shifter 
(WLS) mounted on the rear end of the crystals and ful- 
filled two purposes. The WLS collect the light from the 
crystals over their full cross sectional area and transform 
the wavelength of the detected light to the sensitive range 
of the photodiode. Csl crystals emit light at wavelengths 
between 450 and 610 nm, with a maximum at 550 nm. The 
absorption profile of the WLS exhibits a maximum be- 
tween 520 and 590 nm matching the range of emission of 
Csl. The light emitted by the WLS has a wavelength of 
600 to 700 nm. For this wavelength, the WLS is practically 
transparent and the photodiode is most sensitive. 



Preamplifiers are attached to the photodiodes on the 
crystal module to reduce the background noise. The over- 
all gain stability of the Csl crystals is monitored by a light 
pulser system. For this purpose each WLS is connected via 
an optical fiber to a common Xenon flashlamp operating 
at a repetition rate of 5 Hz. The light is read out by the 
normal calorimeter electronics. Thus, crystal modules not 
working properly can easily be identified. In addition, the 
light pulser system is used to calibrate the high and the 
low range of the 12bit-dual range ADCs. 

The Crystal Barrel was already used successfully for 
seven years at CERN 11 (LEAR 12 ) before it was brought 
to Bonn. Changes of the photon yield in the crystals due 
to radiation damage were not observed. A detailed de- 
scription of the CERN detector can be found in [28]. 

The calibration of the Crystal Barrel was carried out 
after data taking. For each crystal, the ir° peak in the 
two-photon invariant mass was normalized iteratively to 
the nominal mass of m n o = 134.98 McV by taking the 
invariant masses of each pair of two photons. The mass 
resolution is a = 8MeV at m„o and a = 15MeV at m v in 
their 2-f decays, and a = 20MeV at m u (in its 7r°7 decay). 



3.5 Inner Scintillating-Fiber Detector 

The three-layer scintillating fiber detector identifies charged 
particles leaving the target and determines their intersec- 
tion point with the detector [29] . The 40 cm long detector 
is mounted on a 1.8 mm thick aluminum support struc- 
ture. The fibers are 2 mm in diameter; they are positioned 
at mean radii of 5.81cm, 6.17 cm and 6.45 cm, respec- 
tively. The innermost layer corresponds to a solid angle 
of 92.6% of 4-7T, thus covering almost the full solid an- 
gle of the barrel detector. The fibers of the outer layer 
were installed parallel to the z-axis, the fibers of the in- 
ner layer are bent clockwise, the remaining fibers are bent 
anti-clockwise, forming respectively, an angle of —24.5° or 
25.7° with the z-axis. Between each pair of layers, carbon 
cylinders hold the fibers in place. In total, the detector 
consists of 513 scintillating fibers read out via 16-channel- 
photomultipliers. Each scintillating fiber is coupled indi- 
vidually to an optical fiber guiding the signal through 
the backward opening of the barrel to the photomulti- 
pliers. The 3 layers had efficiencies of 94.8 ±0.9% (inner), 
92.9 ± 0.9% (middle), and 88.1 ± 0.8% (outer), respec- 
tively. These values include the geometrical acceptance. 
The probability of two out of three layers having fired was 
98.4 ±0.2% [29]. 



3.6 Trigger and Data Acquisition System 

The first-level trigger demanded a coincidence of a hit 
in one of the tagger scintillators (related to an incoming 
photon) and of hits in two or three different layers of the 
inner scintillating-fiber detector (interpreted as detection 
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of a proton in the final state). The tagger rate was of the 
order of 1-3T0 6 Hz; the coincidence with scintillating-fiber 
hits reduced this rate to 2000 Hz. The second level trigger 
defined the number of contiguous clusters of hit crystals. 
It used a fast cluster encoder (FACE) based on cellular 
logic. The decision time depended strongly on the com- 
plexity of the hit distribution in the Crystal Barrel and 
was typically 4/is. In case of rejection of the event, a fast 
reset was generated, which cleared the readout electronics 
in 5 [is. Otherwise the readout of the full event was initi- 
ated; typical readout times were 5-10 ms. Along with an 
incoming trigger rate of 2000 Hz from the inner detector, 
this led to a dead time of 70% when two clusters were 
required by FACE. For three clusters in FACE, the dead 
time was negligible. 

For the data presented here, events were only recorded 
if they had more than two clusters determined by FACE 
from the pattern of hit crystals. For part of the data the 
minimum number of clusters was set to three. The average 
overall data taking rate of the DAQ was 100 Hz or more, 
depending on whether empty channels were suppressed or 
not. 



4 Data Analysis 

The data presented here were taken from December 2000 
until July 2001 in two run periods with different primary 
electron energies of 1.4 and 3.2 GeV. In the following, we 
refer to these data sets according to their incident electron 
energies. Cross sections were determined separately for the 
two different settings providing a total range of photon 
energies from 0.3 to 3 GeV. This corresponds to a range 
of 7p or p7r° invariant masses from 1.2 to 2.6 GcV/c 2 . The 
7T° is reconstructed from its 7r° — > 77 decay. 



4.1 Photon Reconstruction in the Crystal Barrel 

An electromagnetic shower will, in general, extend over 
several crystals. Such an area of contiguous calorimeter 
modules is called a "cluster" . In order to reduce the con- 
tributions of noise to the cluster energy, only crystals with 
an energy deposit above a threshold of 1 MeV arc consid- 
ered in the cluster finding algorithm. Within a cluster, a 
search was made for local energy maxima with an energy 
deposition above a minimum value of 20 MeV (-Ecls)- A 
crystal containing a local maximum of energy is referred 
to as a central crystal. Any local maximum is interpreted 
as evidence that a (charged or neutral) particle has hit the 
Crystal Barrel. Its energy deposit in the Barrel is called 
PED (particle energy deposit). If a cluster contains only 
one local maximum, the sum over all crystal energies in 
the cluster was then assigned to the energy of the PED 
(.Eped)- The center of gravity x = (0, of the energy 
distribution defines the spatial coordinates of the impact 
point on the Crystal-Barrel surface, and the momentum 



of a photon: 
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where j runs over all crystals in the cluster and the cut- 
off parameter is P — 4.25. A logarithmic weight has a 
stronger impact from small (and more distant) energy con- 
tributions leading to a better reconstruction accuracy of 
the photon direction. The weighting procedure was opti- 
mized performing extensive Monte-Carlo simulations. The 
PED is identified with a photon if it cannot be associated 
in space with the projection of a charged track (see be- 
low) extrapolated to the surface of the Crystal Barrel. 
Otherwise, it is identified as charged particle and here, by 
default, as proton. 

If more than one local energy maximum is found in- 
side a cluster, the energy of PED i is determined from the 
energy deposit of the central crystals and the deposits of 
the up to eight neighbors, summed up to form Eg. The to- 
tal energy of the cluster is then shared between the PEDs 
j = 1, 2, • • • according to the relative magnitude of their 
_Eg-sums: 
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where crystal energies in overlapping regions are split ac- 
cording to the relative energy deposits of the local max- 
ima. 

Monte Carlo simulations showed that the reconstructed 
photon energies differ from the initial values due to energy 
loss in insensitive parts of the detector. A PED-energy- 
and (9-dcpendcnt function derived from MC simulations 
was applied to the reconstructed energies. 

Statistical shower fluctuations can generate additional 
maxima which are called split- off "s and may be misidenti- 
fied as photons. A single photon can create several split- 
offs, so that some good events are not considered in the 
analysis due to a wrong PED multiplicity in the barrel. 
About 5% of all photon showers create such split-offs. 
Split-offs can also be generated by charged particles, e.g. 
by nuclear reactions. The Monte Carlo simulation repro- 
duces split-offs rather precisely. The error originating from 
split-offs is estimated to be less than 2% [30]. 



4.2 Identification of Charged Particles 

The event reconstruction in the Crystal Barrel provides 
direction and energy of a photon (<£•, 0, and E + with er- 
rors). High-energy charged particles like protons and pions 
can traverse the whole length of a crystal module with- 
out depositing all their kinetic energy in the Csl crystals. 
Pions originating in the target center and traversing the 
Crystal Barrel show a minimum ionizing peak at about 
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Fig. 7. Charged particles crossing the inner fiber detector 
produce different hit patterns. A hit in three layers may re- 
sult in one unique intersection point or in three intersection 
points (middle). One inefficiency still yields a defined intersec- 
tion point (right). 

170 MeV. Their kinetic energy cannot be deduced from 
their energy deposit in the Crystal Barrel. The energy de- 
posit of protons is larger and reaches a maximum for those 
stopping at the rear end of the Crystals. Protons stop- 
ping earlier and penetrating protons deposit less energy. 
A given PED energy of a proton corresponds therefore to 
two allowed values for the proton kinetic energy. In this 
analysis, the proton PED energy was not used. 

Identification of charged PEDs was important to per- 
form the analysis even though proton PEDs were ignored 
after identification. Charged particles produce usually single- 
crystal clusters, and the reconstructed angles and <P 
agree in most cases with a crystal center, leading to reg- 
ular spikes in the proton angular distribution. Therefore, 
protons were treated as missing particles in a kinematic 
fit. 

A charged particle traversing the fiber detector can 
fire one or two adjacent fibers in each layer forming a 
British-flag-like pattern which defines the impact point 
(see Fig. 7, left). A single charged particle may also fake 
three intersection points (see Fig. 7, middle). In order not 
to loose these events, up to three intersection points were 
accepted in the reconstruction. One inefficiency (one bro- 
ken or missing fiber) leading to a pattern shown in Fig. 7, 
right, was accepted by the reconstruction routines. 

The accuracy of the reconstruction of the impact point 
was studied using simulations and was determined to ±0.5 
mm in the x- and y- coordinates (representing the resolu- 
tion in <P). The error in the z coordinate was determined 
to be 1.6 mm. Trajectories were formed from the center of 
the hydrogen target to the impact point in the fiber detec- 
tor as well as to all PEDs in the barrel. For all resulting 
combinations, the angle between the trajectories was cal- 
culated. A PED in the Crystal Barrel was identified as a 
charged particle if the angle between a trajectory start- 
ing in the target center and hitting the fiber impact point 
and the trajectory from target center to the Crystal-Barrel 
PED was smaller than as 20° ( = 0.35 rad). This procedure 
is called matching. 

In some cases, there was more than one hit per layer. 
If the hits belonged to adjacent fibers, the clusters were 
assigned to the central fiber (two wires give an effective 
fiber number (2n+l)/2); otherwise, more than one impact 
point was reconstructed. Each impact point was tested 



if it can be matched to a PED in the Crystal Barrel. 
Matched PEDs are identified as protons. Only one pro- 
ton is allowed in the event reconstruction. Events with 
two matched PEDs, possibly due to p7r + 7r~7r° with one 
undetected charged particle, were rejected. 

Protons going backwards in the center-of-mass system 
have rather low momenta. The minimum momentum for a 
proton to traverse the target cell, wrappings, inner detec- 
tor, and Crystal-Barrel support structure and to deposit 
20 MeV in a crystal is 420MeV/c. To reconstruct events 
with protons having smaller momenta, events with two 
PEDs in the Crystal Barrel were taken into account when 
they had a hit in the inner detector not matching one of 
the two PEDs. It was then assumed that the two PEDs 
were photons and that the proton got stuck in the inner 
detector or support structure between scintillation fiber 
detector and the barrel. Thus, proton momenta down to 
260MeV/c were also accessible in the analysis. 



4.3 Event Reconstruction 

In the first step of data reduction, both photons from a 7r 
decay were required to be detected in the Crystal Barrel, 
irrespective of the detection of the proton. Thus, events 
with 2 or 3 PEDs in the Crystal Barrel were selected. If 
a proton candidate was found in the fiber detector, and 
the fiber hit could be matched to a barrel hit, the corre- 
sponding PED was identified as proton and excluded from 
further analysis. The proton was thus treated as missing 
particle. Unmatched PEDs were considered to be photons. 
All events with two photons were kept for further analysis. 
In this way, consistency between the two data sets (with 
and without observed proton) was guaranteed. 

In total, ~45 million events were assigned to these two 
event classes and subsequently subjected to a kinematic 
fit as described in the following. 

Kinematic fits minimize the deviation of measured quan- 
tities from constraints like energy and momentum conser- 
vation. The measured values yt are varied (by 5yt) within 
the estimated errors &i until the constraints arc fulfilled 
exactly. Further constraints are given by the known masses 
of particles like ir° or r\ which are reconstructed from the 
4-momenta of photons. Kinematic fitting improves the ac- 
curacy by returning corrected quantities, which fulfill the 
constraints exactly. A \ 2 value can be calculated from the 
minimization given by 

* - ? (M ■ 

% ^ 7 

In our experiment, the quantities yi are the measured val- 
ues (<l>, 0, v~E). Since a fraction of the energy can be lost 
(but never created) in some material, the distribution of 
E is asymmetric; therefore \f~E~ is chosen as variable which 
exhibits a more Gaussian distribution than E. 

If the yt measurement errors are correlated, then this 
becomes a matrix equation. Each constraint equation is 
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linearized and added, via the Lagrange multiplier tech- 
nique, to the x 2 equation: 



X 



(6y) T V-\6y) + 2a T f(y + <5y) , 



OS) 



where V is the covariance matrix for the measurements 
y = (yi); a is the vector of Lagrange multipliers for the 
constraints f , and Sy denote the required variations. When 
the constraining equations are satisfied and the experi- 
mental error matrix correctly determined, then the differ- 
ence between true values yi and the measurements yi + Siji 
should be of the same magnitude as the errors. If the data 
are distributed according to a Gaussian function, a prob- 
ability distribution for the x 2 with n degrees of freedom 
can be defined by 



P(x 2 ;n) 



2 -n/a 



X 



-2 e -x-/2 



(9) 



which is called the \ 2 probability. To make judgments and 
decisions about the quality of the fit, the relevant quantity 
is the integral 



CL( X 2 ) = 



P(x 2 ;n) d X 2 



(10) 



which is called the confidence level. The x 2 calculation 
depends directly on the errors of the measured values. In 
this case, the errors are the measurement errors of the 
calorimeter hits: the polar angle 0, the azimuthal angle <P 
and the energy E of a calorimeter cluster. If the errors and, 
hence, the covariance matrix are correctly determined, the 
confidence level distribution of the fits should be flat. A 
sharp rise near zero of the confidence level distribution 
can indicate contributions from background events. The 
kinematical fit is thus an ideal method to quantify the 
quality of different final-state hypotheses for an event. 

A pull is a measure of the displacement of the measured 
values to the fitted values normalized to the corresponding 
errors. The measured values yi with errors o~i are corrected 
by shifts 5yi leading to new values yi + Syi which fulfill 
exactly the constraints and have smaller errors a'i. The 
quantity 

Syi 



pull; = 



(11) 



is called pull; it should be Gaussian with unit width and 
centered around zero. Fig. 8 shows the pulls for the Crys- 
tal Barrel quantities (<P, 0, vE) and the photon beam 
energy E 1 . Mean values are all compatible with zero, the 
variances are close to (but slightly above) 1. 

A systematic offset in a measured quantity will move 
the center of its pull distribution away from zero. Using 
this information in addition to kinematic fits treating the 
z vertex as free parameter, an (unwanted) target displace- 
ment of 0.65 cm from the central position towards the tag- 
ging system was detected, in perfect agreement with a 
later position measurement. 

In a preselection, the compatibility of events with the 
hypothesis 

7P - P27 (12) 



was tested in a one-constraint (1C) kinematic fit imposing 
energy and momentum conservation but leaving the pro- 
ton 3-momentum as adjustable quantity. The confidence 
level distribution of the fit is shown in Fig. 9. Above 20%, 
the confidence level distribution is flat in both data and 
Monte Carlo events. At small confidence level CL, the dis- 
tribution increases when CL approaches zero in both dis- 
tributions. In the data, this could indicate the influence 
of background events or of too small errors for a subclass 
of events. In the Monte Carlo simulation shown in Fig. 9, 
only p7T° events were created and no background events. 
Hence the events with small confidence level must be true 
p7r° events. Indeed, inspection of 77 invariant mass distri- 
bution shows that these events are p7T° events with larger 
errors. The assumption made in the kinematic fit that the 
errors do not depend on the kinematics of the event is ob- 
viously wrong. A cut at e.g. 10% confidence level entails 
therefore the risk that events with larger errors, prefer- 
entially events with ir° in forward direction, have only 
a small detection efficiency. This would lead to a loss of 
statistics and, in case of differences between data and sim- 
ulations, to wrong results. 

Here, a cut on the confidence level at 10 -4 was applied. 
This cut rejects all events where the kinematic fit fails 
(CL=0) and very badly measured events. It was checked 
that the rejected events have a 77 mass distribution with 
very few 7r° only. 

Fig. 10 shows 77 invariant mass spectra after the con- 
fidence level cut: (a) for the low-energy run and (b) for the 
high-energy run. There arc approximately 2.44- 10 6 (1.32- 
10 6 ) 7T° events in the low-energy and 572 -10 3 (144- 10 3 ) ir° 





Fig. 8. Deviations between the measured values 4> (a), (b), 
\[E (c) and Ej (d) as the results from a kinematic fit in units 
of the respective measurement errors (1.4 GeV data). After cal- 
ibration, the mean values should be and the variance 1. The 
experimental distributions are compatible with these numbers. 
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events in the high-energy data set. The numbers in paren- 
theses give the number of events with a proton seen by 
the Crystal-Barrel detector. This fraction is larger for the 
high-energy run where the tagged energy range started at 
higher energy, and more protons escaped through the for- 
ward hole. For a fraction of this data, a more restrictive 
trigger excluded events with only 2 detected particles. 

The 7T° meson is observed above an almost negligible 
background at a level of 10 -3 . The remaining background 
under the ir° at high energies was subtracted using side 
bins, which contained typically a few events. Empty-target 
runs were used to determine additional background. Af- 
ter all cuts had been applied, very few events survived in 
empty-target runs. In LH 2 data, only 2 ± 2 % were back- 
ground events not stemming from LH2. 

Cuts in the invariant mass were applied in addition 
to the kinematic fits. The criterion for a 7r° meson was a 
27 mass in the (135±35) MeV/c 2 interval. At high photon 
energies, the cut was widened to the 75 - 175 MeV/c 2 
region. For the low-energy run, a second kinematic fit was 
applied constraining the two-photon invariant mass to the 
7T° mass (two-constraint fit). 



4.4 Monte-Carlo Simulations 

The detector was simulated by CB-GEANT, a Monte Carlo 
program based on the CERN program package GEANT3. 
The geometry of the barrel calorimeter was implemented 
accurately. Only the type-11 crystals (see fig. 6) had to be 
approximated since GEANT3 did not support their non- 
trapezoidal shape. 

A number of small energy correction factors were ap- 
plied to the simulated energy deposits in the Csl crystals. 
Photons of larger energy have a shower profile penetrating 
deeper into the crystals. The efficiency of light collection 
is higher by several per cent at the rear end of a crystal. In 
addition, the light collection efficiency is different for dif- 
ferent crystal shapes. The mean energy loss due to shower 
leakage into support material was determined from the 




(b) 3.2 GeV data 

Fig. 10. Two-photon invariant mass spectra using events 
with two properly reconstructed photons. Note the logarith- 
mic scale. A kinematic fit to the 7p — ► P77 hypothesis was 
performed and a confidence level cut of 10~ 4 applied. The dark 
grey area, defined by an additional mass cut, corresponds to ac- 
cepted events. The low-mass asymmetry of the 7r° peak is due 
to energy overflow in the ADCs for high-energetic photons. It 
is well reproduced by Monte Carlo simulations. 
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Fig. 9. Confidence level distributions (number of events as 
a function of CL) for the one-constraint kinematic fit to the 
hypothesis 7p — > p27 imposing energy and momentum con- 
servation and treating the proton as unseen particle for 1.4- 
GeV data (left) and the corresponding Monte Carlo simula- 
tions (right). 



simulations. These effects required empirical corrections 
of the reconstructed photon energies. 

The fiber detector was implemented as a homogeneous 
cylindrical scintillation counter since simulating fibers with 
helix-shaped bending is not supported by GEANT3. In- 
stead, the fiber number in each layer was calculated from 
the known impact point taking into account the real shape 
of the fibers. Single- fiber inefficiencies as derived from the 
data were taken into account to allow simulations of trig- 
ger efficiencies. The gaps between the cylinders were filled 
with a carbon-based support structure. The plastic-foil 
wrapping was simulated as described in (5.3.1). 

The acceptance was determined from Monte-Carlo sim- 
ulations as ratio of reconstructed and generated events: 



gen,MC 
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Fig. 11. Acceptance as a function of tagger wire number for the 20 cos cm s bins from [—1.0, 0.9], . . . , [0.9, 1.0]. The acceptances 
are plotted at the energy _E 7 of the corresponding wire. The acceptance is calculated using eq. (13). At large energies, the energy 
range covered by one wire is small, and the acceptances show statistical fluctuations due to limited Monte Carlo event numbers. 



The acceptance is shown in Fig. 1 1 for the 20 cos <9 cm bins 
and for both data sets as a function of the wire number. 
Instead of the wires number, the corresponding photon 
energy is plotted on the horizontal axis. A^a_, 2 - t reaches a 
maximum of 60% to 70% depending on the energy range. 
The acceptance is very small for cos(9 cm > 0.9. The cor- 
responding low-energetic protons hardly reach the fiber 
detector and thus the first-level trigger condition is not al- 
ways met. For cos0 cm < —0.9 (for backward 7r°), protons 
simply escape through the forward hole not covered by 
the fiber detector. For high-energy ir° mesons, the opening 
angle between the two decay photons is not large enough, 
and only a single PED is detected. Then the event is lost in 
the trigger and the detection efficiency is smaller at higher 
energies. To avoid biased results due to imperfections of 
the apparatus and/or the simulations, we required the de- 
tection efficiency to exceed a minimum of 5 %. For this 
reason, we rejected the last forward and backward data 
points. 



5 Determination of Cross Sections 
5.1 Basic Definitions 

The unpolarized differential cross sections da/dfl can be 
calculated from the number of data events identified in 



the respective channel using 



dcr 



-2- 



1 



1 



Aotal 



A. 



N lPt An /> 



*2~> 



(14) 



where the quantities are: 



N 7 
Pi 

An 



*2j 



Number of events in (E 7 , cos (9 cm ) bin, 
Acceptance in (i? 7 , cos (9 cm ) bin, 
Number of primary photons in E bin, 
Target area density, 

Solid-angle interval An = 27rZ\(cos @ cm ), 
decay branching ratio. 



The number of events in a (i? 7 , cos (9 cm ) bin comprises 
events with two or three PEDs in which the proton was 
either undetected or detected in the Crystal Barrel. Sum- 
mation of both contributions reduces the necessity of re- 
producing exactly the threshold behavior of low-energetic 
protons in the Crystal-Barrel detector. For part of the 
3.2-GcV data run, three or more hits in the Crystal Bar- 
rel were already required by the hardware trigger (FACE). 
The most forward data point has thus a smaller number 
of events and a larger error. 

The target area density is calculated from the density 
of liquid hydrogen p(LH2) = 0.0708 g/cm 3 and its molar 
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mass M mo i(LH 2 ) = 2.01588 g/mol to be 



5.3 Systematic Uncertainties 



Pt 



2p(LK 2 )N A L 
M mol (LH 2 ) 



2.231 • KT 7 /fjb, (15) 



where N\ = 6.022 • 10 23 /mol is the Avogadro's number, 
I = 52.75 mm is the length of the active target cell, and 
the factor 2 accounts for the two atoms in LH2. 

The 7T° was identified via its decay into 27 which has 
a relative branching ratio of 98.798 %. 

The solid-angle interval is 



AQ = 2?rZ\(cos 6> c 



(16) 



where Z\(cos (9 cm ) = 0.1 gives the bin width of the angular 
distributions, subdividing cos <9 cm into 20 bins. Photon- 
energy bins of about 25MeV width were chosen for the 
1.4 GeV data. The 3.2-GeV data are presented in bins 
of about 50MeV, 100 MeV, and 200 MeV in the intervals 
E 1 e [750, 2300], [2300, 2600], [2600, 3000], respectively. 



5.2 Normalization 

At E e , elsa = 1.4 GeV, the SAID-SM02 model descrip- 
tion of 7p — > p7T° was used for the absolute normalization 
of our measured angular distributions. We minimized the 
function 



ra-l 

E 

i=0 



(x-Ni- Stf 



(S 



(17) 



for each energy channel of the tagging system, where i 
denotes a bin in cos6> cm , iVj are the acceptance-corrected 
number of data events, 5^ are the errors in iVj, and Si 
are the values of the SAID model integrated over the cor- 
responding energy bins. Furthermore, x stands for a nor- 
malization factor including all normalization constants of 
Eq. (14), it is used as free fit parameter. This method 
works well for the energy region between threshold and 
E 1 w 2 GeV, where the SAID model is reliable. The ex- 
tracted photon flux leads also in the channel 7p — > pr] 
to cross sections which are consistent with previous data. 
We estimate the error of the normalization to be ±5%. 
The error given by the statistical spread of all data on 7r° 
photoproduction used by SAID is not taken into account. 

The method of normalizing to a known cross section 
was not applicable to our data at E ei elsa = 3.2 GeV since 
our data extend the currently available world database 
substantially. Instead, we relied on the measured number 
of electrons in the tagger, which is proportional to the pho- 
ton flux. A normalization factor was determined from a fit 
to SAID-SM02 in the energy range up to 1.7 GeV and ap- 
plied to the full energy range. Some systematic deviations 
in the absolute cross section of the order of 10-15% occur 
in the lowest bins, for 800 - 950 MeV. The deviations are 
not caused by uncertainties in detection efficiency. This 
follows from the perfect agreement of the 27 and 37r° de- 
cay channels of the r/ in 7p — > p?y [7]. Thus, we believe that 
within a systematic error of 15 % we can trust this normal- 
ization, even above incoming photon energies of 2 GeV. 



Systematic errors were studied in Monte-Carlo simula- 
tions. The errors, with exception of the normalization er- 
ror, were added quadratically for each point in the dif- 
ferential cross section. The total systematic error is then 
added quadratically to the statistical error. The system- 
atic errors discussed in the following subsection arc in- 
cluded in the shown error bars. The statistical errors are 
small, except for the most forward and most backward 
points in the angular distributions. 

The following effects contribute to the systematic un- 
certainty of the measurements: 



5.3.1 Reconstruction Efficiency 

The reconstruction of neutral mesons and the identifica- 
tion of final states required a sequence of cuts including 
those on the results of kinematic fitting. An overall un- 
certainty of ±5.7% was assigned to the reconstruction ef- 
ficiency as determined in [30]. This error includes uncer- 
tainties due to split-offs. 

5.3.2 Target Position 

The position of the target cell was determined by compar- 
ing results from kinematic fitting (off-zero displacement of 
pull distributions) to Monte-Carlo simulations. The target 
was found to be shifted upstream by 0.65 cm, i.e. into the 
direction of the tagger with respect to the center of the 
Crystal Barrel. The acceptances of the 7p — > p7T° reaction 
were rc-detcrmincd for target shifts of ±3 mm compared 
to the real target position. The variations in the differen- 
tial cross sections that resulted from these changes depend 
on energy and angle are as large as 5% for forward protons 
and ±1% on average. 

5.3.3 Position of Beam Axis 

The position of the photon beam measured by a beam 
profile monitor showed slight variations in time and even 
within an extraction cycle. Analyzing data from the inner 
detector, this shift was found to be less than 3 mm off axis 
at the target position. The acceptances were re-calculated 
assuming a shifted vertex distribution and the resulting 
cross section changes were included as systematic uncer- 
tainty. The errors due to beam shifts are, depending on 
energy and angle, 2% or smaller. 

5.3.4 Material between Target and Crystal Barrel 

All material between target and sensitive detector com- 
ponents was simulated carefully. However, black tape and 
foil was used to wrap the inner detector to get it light- 
tight; its thickness was not exactly known. Small changes 
in the thickness of this material had large effects on the 
resulting detection efficiency for low energy protons and 
on the trigger efficiency (which required at least two layers 
of the scintillation fiber detector to have fired). Assuming 
the worst case of 1 mm of additional or missing material, 
simulations were performed, new acceptances calculated 
and systematic uncertainties in the cross sections deter- 
mined. This error was negligible. 
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5.3.5 Solid Angle 

Events belonging to one bin in cos <9 cm may be recon- 
structed in an adjacent bin. This migration effect was 
studied by unfolding the resolution in cos(9 cm by Monte 
Carlo methods and found not to contribute significantly 
to the results. 

5.3.6 Target Thickness 



calculating mean values for 0.8 < E~ t < 1.3 GeV since the 
final errors arc dominated by common systematic errors. 
For this energy region, we used the 1.4-GeV data because 
of the finer binning and since the tagger worked more re- 
liably in the high energy part where the intensity is lower. 

6.3 Differential Cross Sections dcr/dfl for 7p — *■ p7r° , 
Combined Data Set 



The target thickness does not contribute to the error since 
our cross sections are normalized to SAID. 

5.3.7 Photon Flux 

Our angular distributions are normalized to SAID, as out- 
lined in section 5.2. We assign an error of ±5% (±15%) 
to the normalization of the 1.4 GeV (3.2 GeV) data. These 
errors arc not included in Figs. 12-14. 



6 Experimental Results 

Differential cross sections were calculated separately for 
the two different ELSA energies of 1.4 GeV and 3.2 GeV. 
The data sets published in [6] contained the complete 1.4- 
GeV data set covering the photon energy range from 0.8 
to 1.3 GeV and the 1.3 to 3.0 GeV range from the 3.2-GeV 
data set. 



6.1 Differential Cross Sections dcr/dfl for 

7P — > p7T° at an Electron-Beam Energy of 1.4 GeV 

The differential cross sections are shown in Fig. 12 to- 
gether with the SAID-SM02 model curve. Statistical and 
systematic errors are added quadratically. In general, the 
agreement is impressive. Note, however, that the data are 
normalized to SAID in each energy bin. At the lowest ener- 
gies, we encounter some deviations from SAID. We cannot 
exclude systematic effects beyond those listed in section 
5.3. They may be due to the rather small kinetic energies 
of the recoil protons from the reaction 7p — > p7r° , just 
at the threshold to trigger the inner detector. The high 
tagger rates in this energy region may also be responsible 
for the problems. 



6.2 Differential Cross Sections dcr/df2 for 
7P — > p7T° at 3.2 GeV 



Fig. 13 shows the angular distributions for the reaction 
7P — ► p7T° for the high-energy data set. SAID-SM02 
results are shown for comparison. There is good overall 
agreement between data and SAID. At the smallest ener- 
gies, the agreement is somewhat worse. For this data, one 
normalization factor is used. 

The two data sets are combined by using the 1.4-GeV 
data set up to the maximum possible photon energy of 
1.3 GeV and using the 3.2 GeV data set to cover the range 
from 1.3 GeV to 3 GeV photon energy. We restrained from 



The cross sections from threshold to 3 GeV are shown in 
Fig. 14. Data points are shown as open circles with er- 
ror bars. The solid line represents a fit described below. 
The MAID results [25] for photon energies up to 1 GeV 
are given by dotted lines, results of the new SAID model 
SM05 for photon energies up to 2 GeV as dashed lines. 
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Fig. 12. Differential cross sections for 7p — » p7r° from the 
low-energy data. Statistical and systematic errors are added 
quadratically. The photon flux was determined by normalizing 
the distributions to the SAID-SM02 model. The photon energy 
range (in MeV) is given in each subfigure. The symbols are: ■: 
CB-ELSA, — : SAID-SM02. 
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The agreement between data and the models is equally 
good for both data sets. The angular distributions show 
rapid variations depending on both energy and emission 
angle of the meson, reflecting the contributions from many 
partial waves. 
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Fig. 13. Differential cross sections for jp — ► p7r° from the 
high-energy data. Statistical and systematic errors are added 
quadratically. The absolute normalization was determined by 
fitting the angular distributions in the range from 800 MeV to 
1.7 GeV to the SAID-SM02 model. The photon energy range 
(in MeV) is given in each subfigure. The SAID-SM02 model 
result is shown as solid line. 
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Fig. 14. Previous page: Differential cross sections for 7p — > 
p7T° ; combined data set. Statistical and systematic errors are 
added quadratically. The solid line represents a fit to the partial 
waves described below, the dashed line is the SAID SM05, the 
dotted line the MAID model. 



6.4 Partial Wave Analysis 

The strong variations in the differential cross sections in- 
dicate strong contributions of various resonances with dif- 
ferent quantum numbers. To determine quantum numbers 
and properties of contributing resonances, the differential 
7P — > p7r° cross sections were used in a partial wave anal- 
ysis. The analysis was based on an isobar model. In the s- 
channel, N* and A* contribute to the p7r° final state. The 
background in this channel was described by reggeized t- 
channel p (u) exchange and by baryon exchange in the it- 
channel. In addition s-channel Born-terms were included 
in the fits. Details on the partial wave analysis can be 
found in [9,10]. The CB-ELSA data on 7 p -»■ P 7r° and 
7P — * prj [6,7] were included as well as additional data 
sets from other experiments: Mainz-TAPS data [31] on 
rj photoproduction, beam-asymmetry measurements of ir° 
and r\ [23,32,33], and data on 7p — > nir + [34]. The high 
precision data from GRAAL [23] do not cover the low 
mass region; therefore we extracted further data from the 
compilation of the SAID database [32]. Data on photo- 
production of K+A, K+E°, and K°E+ from SAPHIR [35, 
36] and CLAS [37], and beam asymmetry data for K + A, 
K + E° from LEPS [38] were also included in the analy- 
sis. This partial wave analysis is much better constrained 
than an analysis using the 7p — > p7r° data only. The main 
results on baryon resonances coupling to p7r° and p?7 are 
discussed in [10], those to K+A and K + E° are documented 
in [11]. 

To describe the different data sets, 14 N* resonances 
coupling to N7r, Ntj, KA, and KE and 7 A* resonances 
coupling to N7r and KE were needed. Not all included 
resonances contribute to the p7T° final state. Most reso- 
nances were described by relativistic Breit-Wigner ampli- 
tudes. For the two Sn resonances at 1535 and 1650 MeV, 
a four-channel if-matrix (N-7T, N77, KA, KE) was used. 

The total cross section and the main results of the 
PWA are briefly discussed in the next section. 



6.5 Total Cross Section and Results of the Partial 
Wave Analysis 
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Fig. 15. Total cross section (logarithmic scale) for the reaction 
7P — > p7r obtained by integration of angular distributions of 
the CB-ELSA data and extrapolation into forward and back- 
ward regions using our PWA result. The solid line represents 
the result of the PWA. Four individual contributions to the 
cross section are also shown. 



The first resonance region is the dominant structure of 
Fig. 15. It is due to excitation of the A(1232)P33. There is 
strong destructive interference between A(1232) P33, the 
P33 nonresonant amplitude, and u-channel exchange. The 
N(1440)Pn Roper resonance provides only a small contri- 
bution of about 1-3% compared to the A(1232)P33. In the 
second resonance region, the N(1520)Di3 and N(1535)Sn 
resonances yield contributions which are shown as thin 
lines in Fig. 15. 

The third bump in the total cross section is due to 
three major contributions: the A(1700)D33 resonance pro- 
vides the largest fraction (~35%) of the peak, followed by 
N(1680)Fi 5 (~ 25%) and N(1650)S U (~ 20%). In addi- 
tion the A(1620)S 3 i (-7%) and N(1720)Pi 3 (~6%) res- 
onances arc required. In the fourth resonance region, the 
A(1950)F37 contributes ~ 41% to the enhancement and 
A(1920)P 33 is identified with ~35%. Additionally, the fit 
requires the presence of A(1905)F 35 and A(1940)D 33 . The 
high-energy region is dominated by p(to) exchange in the 
f-channel as can be seen by the forward peaking in the 
differential cross sections. 



From the differential cross sections, the total cross section 
was determined by integration. The integration was per- 
formed by summing over the differential cross sections of 
Fig. 14 and using extrapolated values from the fit for bins 
with no data. In the total cross section, shown in Fig. 15, 
clear peaks are observed for the first, second, and third 
resonance region. The fourth resonance region exhibits a 
broad enhancement at W about 1900 MeV. The decom- 
position of the peaks into partial waves and their physical 
significance will be discussed below. 



6.6 Differential Cross Sections dcr/dt 

The partial wave analysis assigns a fraction of the total 
cross section to i-channel exchange which increases with 
energy. This was already expected from the strong rise of 
the differential cross sections towards cos(0 cm ) ~ 1 in the 
higher energy bins. 

Fig. 16 shows differential cross sections as functions of 
the squared four-momentum transfer between the initial 
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Fig. 16. Differential cross sections da/dt for 7p — t pir° . In 
the range of small values of \t — t m i n \ the data were fitted with 
an exponential function e a+b ' t_tmi "' distribution. 

photon and the ir° in the final state. Even though these 
plots do not provide new information compared to dc/dil, 
they are shown here to emphasize the exponential fall-off 
of the cross section at low t. The differential cross sec- 
tions are plotted against \t — t min \, the difference between 
the actual momentum transfer and the minimal momen- 
tum transfer imposed by kinematics. By definition, the 
four-momentum transfer is always negative and can range 
between £ m i n and t max , which for 7p experiments is given 

by 



^min/max 



(S + TO 2 - TO 2 ,) 2 

is 



(18) 

2 



Note that a = mz + 2m p E~ / and that i m i n corresponds to 
forward, i max to backward production of the ir° . 

The squared momentum transfer is related to the emis- 
sion angle of the pion in the center-of-mass system by 



'■\J(s- (m„o +m p ) 2 ) 



(s — (m^o — trip) 2 ) sin 



2 ^cms 



(19) 



Differential cross sections da/dfi are related to da/dt by 
da da dt 



Using (20), the relation of the two different differential 
cross sections is given by: 



da 

dF 



Air s da/dU 



(s - rri^)yj{s- (m p + m 7r o) 2 )(s - (m p - m^o) 2 ) 



(21) 

For small four-momentum transfers, an exponential 
function is fitted to the distributions shown in Fig. 16. 
The cross sections fall off according to 



der 
dt 



oc exp(a + b ■ \t — t B 



(22) 



where the slope parameter b has a negative value. 

For large t, the cross sections do not exhibit the be- 
havior that would be expected if we had only t-channel 
exchange of mesons. There are structures due to other 
phenomena, in particular due to formation of s-channel 
resonances. The range, in which the distributions can be 
fitted by an exponential is not well defined. Nevertheless, 
we show in Fig. 17 the slope parameter — b of Eq. 22 from 
the fits as a function of incident photon energy E~ t . With 
increasing energy, the slope rises to a maximum value at 
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Fig. 17. Slope parameter for 7p — ► p7r° from a fit using 
e a+b|t-t min | ^ ^ e differential cross sections da/dt. 



E 1 ~1. 5-1. 6 GeV and then decreases again. At large en- 
ergies, above 2.2 GeV, the determination of the slope pa- 
rameter becomes somewhat arbitrary, since the most for- 
ward data points at cos cm = 0.85 are missing. For low 
energies, the slope parameter cannot be determined from 
the data and can even adopt positive values. Obviously, 
resonance production is dominant at these energies. The 
results should hence be interpreted with care. In particu- 
lar, the turnover at E 1 < 1.5 GeV is deduced from very 
few data points. 



7 Summary 

We have reported a measurement of unpolarized differen- 
tial cross sections of the reaction 7p — > p7r° in the photon 
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energy range from 0.3 GcV to 3.0 GeV, thus completely 
covering the baryon resonance region. Above 0.8 GeV, the 
data supersede previous experiments both in solid angle 
coverage and in statistics. In the mass range above the 
first resonance region, strong variations of the differential 
cross sections as functions of photon energy and decay 
angle in the rest frame of the intermediate state suggest 
that many partial waves contribute to 7r° photoproduction 
off protons. This is confirmed by a partial wave analysis 
which, beyond non-resonant contributions, identifies a sc- 
ries of resonances contributing to the 7p — > p7r° reaction 
channel. At large photon energies, the angular distribu- 
tions show a forward rise in 7r° direction signaling dom- 
inance of i-channel exchange contributions at high ener- 
gies. The data provide valuable input to coupled-channel 
isobar-model descriptions of baryons and their excitations. 
A partial wave analysis reveales a rich spectrum of reso- 
nant and non-resonant contributions. 
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